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ABSTRACT 

We study the galactic-scale triggering of star formation. We find that the 
largest mass-scale not stabilized by rotation, a well defined quantity in a rotat- 
ing system and with clear dynamical meaning, strongly correlates with the star 
formation rate in a wide range of galaxies. We find that this relation can be un- 
derstood in terms of self-regulation towards marginal Toomre stability and the 
amount of turbulence allowed to sustain the system in this self-regulated quasi- 
stationary state. We test such an interpretation by computing the predicted 
star formation rates for a galactic interstellar medium characterized by lognor- 
mal probability distribution function and find good agreement with the observed 
relation. 

Subject headings: galaxies: disk instabilities - galaxies: formation - star forma- 
tion: general 



1. Introduction 

Observations of normal spiral galaxies by Schmidt (1959) first suggested that their star 
formation rates (SFRs) scale with their global properties. This relation was extended to 
other galaxies with higher SFRs, such as the nuclear regions of spiral galaxies and Ultra 
Luminous InfraRed Galaxies (ULIRGs) by Kennicutt (1998). These observations have lead 
to an empirical law for star formation, often called the Kennicutt-Schmidt (KS) Law: 

Sstar OC E^^^^ , (1) 

where Sgas and Sgtar are the gas surface density and SFR per unit area, respectively. 

Since star formation is a local process that happens on subparsec scales, the correlation 
with galactic scale (> 1 kpc) quantities, such as averaged Sgas, suggests the existence of a 
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physical connection between galactic and subparsec scales. Motivated by the observed KS 
law, several authors have tried to find a scenario in which a global/large scale property of 
galaxies could trigger and/or regulate star formation (Quirk 1972; Wyse 1986; Kennicutt 
1989; Elmegreen 2002; Li et al. 2005). The KS law may also be explained in terms of 
processes that are primarily local, within star-forming clouds (e.g. Krumholtz et al. 2009 
and references therein). 

From the study of gravitational instabilities in disks, the 'turbulent' Toomre parameter 
Qtm-b = v™J.'^/t/7rGSgas (or related parameters such as the star formation threshold Sent; 
Kennicutt 1989) arises as a natural candidate for a key triggering parameter. However, the 
average Qturb in a galaxy is observed to be close to 1, in galaxies ranging from local spirals 
(Martin & Kennicutt 2001) to starbursts such as ULIRGs (Downes & Solomon 1998). Since 
observed values of Qturb (or Sgas/Eprit) range over at most a factor of a few, the observed 
range in SFRs per unit area of seven orders of magnitude is difficult to explain solely in terms 
of this threshold. Therefore, further investigation is required to determine what controls the 
SFRs on galactic scales. 

An important point, not generally addressed, is that, a disk at the condition of marginal 
Toomre stability can have a range of possible self-regulated states. For example, the nuclear 
disk in a starburst with Q ~ 1 is much more turbulent than the disk of a normal spiral 
galaxy with Q ~ 1. The goal of this work is to study which galactic property triggers this 
more turbulent behavior for some galaxies and why their SFRs can be orders of magnitude 
higher than in more 'quiescent' spiral galaxies. 

This work is organized as follows. We first review gravitational instability analysis in 
order to introduce the largest scale not stabilized by rotation, followed in §2 by a discussion 
of the correlation found between this largest scale and the SFR. Section 3 presents a physical 
interpretation of the correlation found, in terms of self-regulation due to feedback processes. 
In §4 we test such an interpretation by comparing predictions against the observed SFRs. 
Finally in §5, we summarize the results of this work. 



2. The Maximum Scale not Stabilized by Rotation 

In order to introduce the largest unstable length scale in galactic disks, we first review 
some standard results from gravitational instability analysis (Toomre 1964; Goldreich & 
Lynden Bell 1965). For one of the simplest cases of a differentially rotating thin sheet or 
disk, linear stability analysis yields the dispersion relation for small perturbations (Binney 
& Tremaine 2008) uj"^ — — 27rGEgas|A;| -\- k'^C^, where Cg = •>/ dP/dE is the sound speed. 
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E is the surface density, and k is the epicychc frequency given by k^{R) = R%- + 4^^ 
(n being the angular frequency). The system becomes unstable when < 0, which is 
equivalent to the condition Q < 1, where Q is the Toomre parameter and is defined as 
CsK/vrGSgas • In such a case there is a range of unstable length scales limited on small 
scales by thermal pressure (at the Jeans length Ajcans = Cg/^^gas) and on large scales by 
rotation (at the critical length set by rotation, Arot = 47r^GSgas/fi;^). AH intermediate length 
scales are unstable, and the most rapidly growing mode has a wavelength 2 Ajeans (Binney & 
Tremaine 2008; Escala & Larson 2008). 

The maximum unstable length scale in a disk, Arot, is a robust quantity because it 
depends only on the surface density and epicyclic frequency of the disk and not on smaller 
scale physics. Such a length scale has an associated characteristic mass, defined as equal to 
Epas(Arot/2)^, which can be expressed as: 

Mrot = -—T^ ■ (2) 



On the other hand, due to the complex structure and dynamics of the real interstellar 
medium (ISM) in galaxies, which cannot be described by a simple equation of state, there is 
not a well-defined Jeans length at intermediate scales. Therefore, there is no real lower limit 
on the sizes of the self-gravitating structures that can form until the thermal Jeans scale is 
reached in molecular cloud cores (Escala & Larson 2008). 

The combination of the observed correlations of SFRs with galactic properties, and the 
fact that this largest scale not stabilized by rotation is the only well defined galactic scale in 
the gravitational instability problem, is our motivation for exploring a possible link between 
SFR and this characteristic galactic scale. 



2.1. The Law 

In order to test the existence of a link between the largest scale not stabilized by rotation 
and the SFR in galaxies, we will check whether the mass-scale defined by rotation (Eq. 2) 
correlates with the SFR. For a rotationally supported system, the average of this mass-scale 
can be expressed in terms of quantities such as the gas mass and gas fraction (Escala & 
Larson 2008): 

U....3.1,^M,^{^y, (3) 

where Mgas is the gas mass in the disk and rj — Mgas/Mdyn is the ratio of the gas mass to 
the total enclosed dynamical mass within the gas radius (this varies from the disk radius for 
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spiral galaxies to the radius of the nuclear starburst disk/ring in ULIRGs). This expression 
(Eq. |3]) has the advantage that it reduces the scatter due to error propagation compared to 
the original formulation (Eq. 

In Figure [1] we plot the maximum mass scale not stabilized by rotation estimated from 
Eq. [3] against the measured SFR (per unit area) in those galaxies. We plot normal spirals 
as star symbols, the nuclear gas in normal spirals as filled circles, and ULIRGs as open 
circles. For the computation of Mj-ot we have considered only molecular gas masses and 
molecular gas fractions. This is because molecular gas should be intimately related to the 
SFR because it is this gas which eventually forms the stars in Giant Molecular Clouds 
(GMCs). Error bars displayed in Figure [1] show uncertainties estimated by error propagation 
from the uncertainties found in the literature. 

The information for each data point plotted in Figure [T] is listed in Table 1, together 
with a list of references to the works in the literature from which the values were taken. For 
the spiral galaxies, Sj, is estimated from the H„ luminosity, the gas masses are estimated 
by the CO luminosity, and the dynamical masses are estimated using the method listed in 
Table 1. For the nuclear gas in normal spirals, the gas masses are estimated using the CO 
luminosity, the dynamical masses from rotation curves and is estimated using method 
listed in Table 1. Finally, for ULIRGs is estimated from the far-infrared luminosity, gas 
masses are estimated from the CO luminosity and dynamical masses are estimated using 
rotation curves. 

Figure [1] shows a clear correlation between Mrot and the SFR per unit area, which 
supports the idea that this threshold mass has a relevant role in the triggering of star 
formation on galactic scales. The solid black line in Figure [1] shows a least-squares fit to the 
points in the figure and corresponds to a star formation law of S^, oc M^^^, with a scatter of 
0.21 dex. This relation has a level of scatter comparable to the typical scatter found for the 
KS law. Moreover, this is the only correlation of the SFR with a galactic quantity with a 
clear dynamical meaning in terms of stability analysis, and therefore with a clear role in the 
star formation problem. 

In summary, the correlation between the SFR and the maximum unstable mass defined 
by rotation is indeed observed in galaxies, over a range that spans almost eight orders of 
magnitude in SFR per unit area. 
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Fig. 1. — The star formation rate plotted against the critical mass-scale defined by rotation 
Mfot, estimated from Eq. |3] using measured quantities in those galaxies. The open circles 
show data for nuclear starburst disks, the star symbols show for normal spirals galaxies and 
the filled circles show for nuclear gas in spirals. The solid line corresponds to oc M^^'^. 
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Table 1: Galaxy Physical Parameters 



Galaxy 


s* 


Mh2 


Mdyn 


Mh, 


References 






IO^Mq 


10^°Mq 




NGC6946 


0.0132 


3.3 


19"^ 


0.02 


(1), (3) 


NGC4419 


0.00404 


4 


9*- 


0.044 


(2), (4), (5) 


NGC4535 


0.00017 


6.918 


26.3'= 


0.0263 


(2), (4), (6) 

V'^/' V^/' V'-'/ 


NGC5033 


0.00229 


6.708 


37.1'' 


0.018 


(1), (7), (8) 


NGC4254 


0.01995 


13.427 


51.6'' 


0.019 


(1), (4), (9) 


NGC7331 


0.004677 


4.94 


21.5'' 


0.023 


(1) (8) (10) 
y-^ji V"/) V-^"^/ 


NGC4303 


0.018197 


10 


25" 


0.04 


(1), (4), (11) 

v-^/i v^/i y-^-^j 


NGC4647 


0.0069 


2.63 


14.6= 


0.018 


(1), (4), (6) 


NGC4654 


0.0087 


3.236 


20.88= 


0.0155 


(1), (4), (6) 


NGC4321 


0.00851 


14.791 


64.3= 


0.023 


(1)> (4), (6) 


NGC4501 


0.006166 


9.77 


69.3= 


0.0141 


(1), (4), (6) 


NGC4689 


0.00417 


3.09 


13.8= 


0.0224 


(1), (4). (6) 

V /' V /' y^/ 


NGC253 


17.378'' 


0.14 


0.04 


0.35 


(1), (12) 


NGC1614 


61.6595'* 


2 


1.3 


0.154 


(1), (13) 


NGC470 


0.7892= 


0.2 


0.067 


0.299 


(14) 


NGC4102 * 


9.06539= 


0.8 


0.3 


0.26 


(14) 


NGC4102 ** 


3.00654" 


1.4 


1.4 


0.1 


(14) 


NGC3504 * 


1.50835" 


0.8 


0.4 


0.2 


(14) 


NGC3504 ** 


2.8349= 


1.2 


1.2 


0.1 


(14) 


NGC4536 


0.5175= 


1.2 


0.8 


0.15 


(U) 


NGC3351 


0.4421^ 


0.2 


0.117 


0.171 


(14) 


NGC3627 


0.1698'' 


2 


6 


0.03 


(1). (15). (16) 


NGC6240 


74.131'' 


3.1 


1.55 


0.2 


(1), (17) 


NGC5005 


0.31079 


1.03 


2 


0.05 


(2), (18) 


IRAS00057 


95.1 


1.4 


1.4 


0.1 


(19) 


IRAS02483 


63.6 


1.9 


2.235 


0.085 


(19) 


IRAS10565 


104.1 


4 


1.739 


0.23 


(19) 


Mrk231 


492 


1.8 


1.286 


0.14 


(19) 


Arpl93 


487.1 


2.6 


1.625 


0.16 


(19) 


Arp220 disk 


166.5 


3 


1.875 


0.16 


(19) 


Arp220 total 


285.4 


5.2 


3.467 


0.15 


(19) 


Arp220 west 


3555.1 


0.6 


0.15 


0.4 


(19) 


Arp220 east 


905.7 


1.1 


0.314 


0.35 


(19) 


IRAS17208 


413.4 


6.1 


2.346 


0.26 


(19) 


IRAS23365 


164.8 


3.8 


3.455 


0.11 


(19), (20) 



* @300pc ** @1300pc 

a) HI rot. curves b) CO rot. curves c) M^y^ — Lh relation d) Lpm e) Lrc f) LBr7 g) Lho 



(1) Kennicutt (1998); (2) Komugi et al. (2005); (3) Crosthwaite & Turner (2007); (4) Young et al. 
(1996); (5) Kenney et al. (1989); (6) Decarli et al. (2007); (7) Perez-Torres and Alberdi (2007); 
(8) Heifer et al. (2002); (9) Sofue et al. (2003); (10) Thilker et al. (2006); (11) Schinncrcr ct 
al. (2002); (12) Maucrsberger et al. (1996); (13) Alonso-Herrero (2001); (14) Jogee, Scoville and 
Kenney (2005); (15) Reuter et al. (1996); (16) Warren et al. (2010); (17) Engel et al. (2010); (18) 
Sakamoto et al. (2000); (19) Downes and Solomon (1998); (20) Murphy et al. (1996) 



-7- 



3. A Physical Interpretation in Terms of Self-Regulation 

It is well established that galactic disks are observed to be globally in equilibrium, for 
galaxies ranging from local spirals (Martin & Kennicutt 2001) to starbursts such as ULIRGs 
(Downes & Solomon 1998), with measured Toomre Q parameters close to 1, in the case 
of a 'turbulent' version of the Toomre parameter: Qturb — v[™b'^/7rGS, where v[™^j^ is the 
observed velocity dispersion in the gas. This state close to stabihty (Qturb ~ 1) was first 
suggested by Goldreich & Lynden Bell (1965) to be due to a self-regulation feedback loop; 
if Qturb >> 1, in the absence of heating driven by instabilities the disk will cool rapidly and 
the system will eventually become unstable, while if Qturb << 1 then instabilities and star 
formation feedback will be so efficient that enough turbulence will be produced to 'heat' the 
disk towards Qturb ~ 1- 

While most galactic disks are close to marginal Toomre stability, some, such as nuclear 
disks in starbursts, can be much more turbulent than the disks of normal spiral galaxies. 
The reason is that although disks are all at Qturb ~ 1- they can have self-regulated states 
with different levels of turbulence (vj;™^). This can be easily visualized by considering that 
the condition Qturb ~ 1 implies a mass scale of Mrot ~ ^['^turbl'^ using Eq 2. Since the 
epicyclic frequency k varies only between Q and 2Q for centrally concentrated disks, in a 
disk supported vertically by turbulence — R/vcirc ~ H/v[™|,) it is straightforward to 
derive a mass scale of 

Mrot - ^mZ.? (4) 

for a median epyciclic frequency of = 3/2 Q. Since the disk scale-height H is a monoton- 
ically increasing function of v[^j, for a disk supported vertically by turbulence, the velocity 
dispersion of turbulent motions follows yI^X °^ with 77 > 0, for a disk with Qturb ~ 1- 
Therefore, since some disks have a larger mass-scale not stabilized by rotation Mrot, their 
large-scale conditions require that feedback processes produce more turbulence in order to 
achieve Qturb ~ 1- 

The existence of disks with self-regulated states of different v^™^ is particularly im- 
portant because it is believed that turbulence has a role in enhancing and controlling star 
formation (Elmegreen 2002; Krumholz & McKee 2005; Wada & Norman 2007). In its sim- 
plest form (proposed by Elmegreen 2002), the SFR depends on the probability distribution 
function (PDF) of the gas density produced by galactic turbulence, which appears to be 
lognormal in simulations of turbulent molecular clouds and the ISM (Ballesteros-Paredes & 
Mac Low 2002; Padoan & Nordlund 2002; Kravtsov 2003; Mac Low et al. 2005; Wada & 
Norman 2007; Wang & Abel 2009). Moreover, numerous numerical studies support the fact 
that the dispersion of the lognormal PDF is determined by the rms Mach number of the 
turbulent motions (Vazquez-Semadeni 1994; Padoan et al. 1997; Padoan & Nordlund 2002; 
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Federrath et al. 2008, 2010). Therefore, it is expected that the SFR in galaxies scales with 
the velocity dispersion of turbulent motions. From this we can infer, at least qualitatively, 
that a galaxy with a higher Mrot should have a higher star formation activity. 

In summary, as the largest mass-scale not stabilized by rotation Mj-ot increases for a given 
disk, its self-regulated state has an increasingly turbulent ISM (v[™j, oc M^ot with rj > 0, for 
Qturb ~ !)• A self-regulated state with higher v™'^, (which itself controls and enhances star 
formation) implies higher star formation activity, and therefore the correlation between the 
critical mass-scale defined by rotation and the SFR is expected. 



4. Galactic SFR for an ISM cheiracterized by a Lognormal PDF 

In order to quantify the arguments given above, in this section we will compute the 
predicted SFR for an ISM dynamically controlled by turbulence in which the rms velocity 
dispersion of turbulent motions vj;™^ is determined by Mj-ot- We will follow an analysis 
analogous to Wada & Norman (2007), starting with the assumption that the density PDF 
of the multiphase ISM in a galactic disk can be represented by a single lognormal function: 



= ^^exp 
\/2na 



Hp/po) 

2ct2 



dlnp , (5) 



where po is the characteristic density scale and a is the dispersion of the lognormal PDF. 
Although there is evidence for deviations from a lognormal function in the tails of the density 
PDF (Scalo et al. 1998; Federrath et al. 2010), for simplicity we neglect any higher order 
correction. 

If the star formation occurs only in regions whose density is higher than a critical value 
{p > Pc), the SFR per unit volume on a global scale is given by: 

p,^e,(Gp,y/X{p)y , (6) 

where Cc is the efhcency of star formation, Sc — pc/ Po is the critical density contrast for star 
formation, (p)y — poe'^^^^ is the volume- average density and fc = 0.5 1 — Erf ^ ^"^J^^ j is 
the mass fraction of gas whose density is higher than pc (for a derivation see §3 of Wada 
& Norman 2007). Numerous numerical studies have claimed that, in addition to the ap- 
proximately lognormal form of the PDF, the dispersion a is determined by the rms Mach 
number A^rms = vJ;™b/Cs, where Cs is the sound speed (Vazqucz-Semadeni 1994; Padoan et 
al. 1997; Padoan & Nordlund 2002; Federrath ct al. 2008, 2010). When a log-normal PDF 
is assumed, the a — A^rms relation can be expressed as (Padoan et al. 1997): 

a'^Hl + h'MLs), (7) 
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where b is a parameter that varies from 0.3 to 1 and depends, for example, on the compressive 
to solenoidal modes of the turbulence forcing (Federrath et al. 2010). Since the forcing mode 
(and therefore b) may vary across regions of the ISM (Federrath et al. 2010), and since the 
sound speed certainly does vary across the ISM, it is convenient for our analysis to define a 
volume-averaged, b- weighted, sound speed Cs = (if )y • 

Using equations HI |6l and El we can write the SFR per unit area (E* = H) as 



1.72 X 10" 



Up 







pc^yr 



tc^c 



Po 



MqPc- 



3/2 



1-Erf (z(4,MrotH'iCg2^ 



X 



X 



1.95; h 



kpc 



and 



z(5c,M,otH-iq 



IO^Mq 



\nS - 21n 1^1 + sGMrotA 



2 ( In ( 1 + 3GMgrt 



1/2 



(9) 



Figure 2 shows the comparison of the predicted correlation between and Mj-ot from Eq 
[8] against the data listed in Table 1, for several values for the model parameters po? (^c, Cs. 
The thick solid lines in Figures 2a, b, c and d represent the predicted for model pa- 
rameters Po = lM0pc~^,ec = 0.01,5c = 10'^, andCs equivalent to a temperature of 250 
K. The other four curves in Figure 2a show variations in the predicted SFRs for po = 
10"^ 10~^ 10^ and 10^ Mopc'^. The four curves in Figure 2b show SFRs for = 1, 10^^ 10"^ 
and 10~^. In Figure 2c the curves show SFRs for 6c = 10, 10^, 10'^ and 10^. Finally, the curves 
in Figure 2d show SFRs for Cs equivalent to temperatures of T=25, 75, 750 and 2500 K. 

From Figure 2 it can be concluded that with a single set of parameters (thick solid lines) 
Eq [8] is able to successfully reproduce the whole correlation, in contrast with the analogous 
work of Wada & Norman (2007), which did not reproduce the K-S law with a single set 
of parameters and relied on a variable star formation efficiency ec (between normal and 
starburst galaxies) in order to reproduce the observed data. 



5. SUMMARY 



In this Letter we have studied the role of the largest mass-scale not stabilized by rotation 
in galactic disks in triggering star formation activity in galaxies. 



10" 10= 10' 10» 10' 10« 10' 10» 10" 



Fig. 2. — Predicted correlation between and Mj-^t from Eq[8]against the data listed in Table 
1, for several values for the model parameters po, ec, Sc, (Cs). The thick solid lines in all figures 
show the model with po = IMopc"^, ec = 0.01, 5c = 10^ and Cs equivalent to a temperature of 
250 K. a) The other lines show the variations in SFRs (per unit area) for log(po/M0pc~^) = 
—2, —1, 0, 1 and 2. b) The lines show SFRs for logec = 0, —1, —2, —3 and —4. c) Curves show 
SFRs for log^c = 1,2, 3, 4 and 5. d) Curves show SFRs for Cs equivalent to temperatures of 
T=25, 75, 250, 750 and 2500 K. 
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We find that a relation between the largest mass-scale not stabilized by rotation and the 
SFR is observed in galaxies ranging from ULIRGs to normal spirals. This relation has a level 
of scatter comparable to the Kennicutt-Schmidt Law, and is the only known correlation of 
the global SFR with a quantity with clear dynamical meaning in terms of stability analysis. 

We give a physical interpretation for the existence of such a correlation in terms of 
self-regulation. In a given disk, as the largest mass-scale not stabilized by rotation increases, 
its self-regulated quasi- stationary state has an increasingly turbulent ISM (vj;™^ oc M[?q^ with 
> 0, for Qturb ~ 1). Therefore, the role of the critical mass-scale in a disk is to define the 
amount of turbulence allowed to be in quasi-stationary equilibrium. Since a sclf-rcgulatcd 
state with higher v[™|, enhances a higher star formation activity, we expect the existence of 
a correlation between the mass-scale for global stability and the SFR. 

We check the validity of this self-regulation scenario by computing the predicted SFR 
for an ISM dynamically controlled by turbulence in which the rms velocity dispersion of 
turbulent motions Vj™|, is determined by Mrot- We find good agreement between the predicted 
and observed SFRs. 

I would like to thank Richard Larson for valuable comments on an early version of the 
draft, Catherine Vlahakis for proofreading this manuscript and the referee. Brant Robert- 
son, for a constructive report. I'm indebted to Fernando Becerra for performing the error 
analysis and graphical display. I also acknowledge partial support from the Center of Excel- 
lence in Astrophysics and Associated Technologies (PFB 06), FONDECYT Iniciacion Grant 
11090216 and from the Comite Mixto ESO-Chile. 
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